Gapped and gapless spin liquid phases on the Kagome lattice 
from chiral three-spin interactions 
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We argue that a relatively simple model containing only SU(2)-invariant chiral three-spin interactions on a 
Kagome lattice of S = 1/2 spins can give rise to both a gapped and a gapless quantum spin liquid. Our 
arguments are rooted in a formulation in terms of network models of edge states and are backed up by a careful 
numerical analysis. For a uniform choice of chirality on the lattice, we realize the Kalmeyer-Laughlin state, 
i.e. a gapped spin liquid which is identified as the v = 1/2 bosonic Laughlin state. For staggered chiralities, 
a gapless spin liquid emerges which exhibits gapless spin excitations along lines in momentum space, a feature 
that we probe by studying quasi-two-dimensional systems of finite width. We thus provide a single, appealingly 
simple spin model (i) for what is probably the simplest realization of the Kalmeyer-Laughlin state to date, as 
well as (ii) for a non-Fermi liquid state with lines of gapless SU(2) spin excitations. 



Quantum spin liquids are elusive phases of matter where 
spins do not freeze into local ordering patterns even at zero 
temperature [ ] . Typically, such spin liquids come in two gen- 
eral flavors, namely gapped topological spin liquids and gap- 
less spin liquids [ ]. Key properties of topological spin liq- 
uids are that they have ground-state degeneracies depending 
on the topology of the system and, when confined to two spa- 
tial dimensions, they host quasiparticle excitations with ex- 
otic anyonic particle exchange statistics. Many examples of 
such systems are known both with [ , ] and without [ ] time- 
reversal symmetry. Systems in the other broad class, gap- 
less spin liquids, are characterized by gapless excitations that 
are not Goldstone modes arising from spontaneous breaking 
of a symmetry. These states can be further divided into two 
groups: One where gapless excitations occur at singular points 
in momentum space - such as, e.g., the case of algebraic spin 
liquids [ ], the celebrated Kitaev honeycomb model [7] and 
variations thereof [8]. In the other group, the gapless modes 
arise on surfaces in momentum space [8-12], generalizing the 
notion of a Fermi surface to systems with bosonic character, 
frequently referred to as Bose surfaces [12]. However, for 
the latter group of states, few results [12-19] firmly establish 
a connection between such states and microscopic models, 
or identify the fundamental theory characterizing the gapless 
modes on the Bose surface of such non-Fermi liquids. 

Here, we explore the phases that emerge when promoting 
the scalar spin chirality 

Xijk = Si ■ {Sj X S k ), (1) 

which is SU(2) invariant, but breaks both parity and time- 
reversal invariance, to an interaction between spins on a lattice 
built from triangles, 

H = ^ JijkXijk Jijk — ±1 (2) 

t,i,feeA,v 

where i, j, k are always ordered clockwise around a trian- 
gle. We show that depending on the assignment of the chi- 
ralities Jij k to the triangles of the lattice, this interaction can 



give rise to either a chiral topological spin liquid or a gapless 
non-Fermi-liquid state with lines of gapless spin excitations 
in momentum space. 

The existence of gapped chiral spin liquid phases, which are 
obtained as generalizations of fractional quantum Hall states 
to spins, was first hypothesized by Kalmeyer and Laughlin [ ] . 
Soon after this proposal, it was suggested that the scalar spin 
chirality ( 1 ) could be relevant for stabilizing such a spin liquid 
phase [20, 21]. However, making a stringent connection be- 
tween a microscopic Hamiltonian and the sought-after chiral 
spin liquid state has remained elusive for many years. Some 
progress has recently been made in constructing somewhat ar- 
tificial spin Hamiltonians [22-25] involving decorated lattices 
or long-range interactions. Furthermore, an alternative per- 
spective arose from topological fiat band models [26-29]. 

Here, we develop a powerful perspective rooted in the 
physics of network models of edge states akin to the Chalker- 
Coddington network model for the integer quantum Hall tran- 
sition [ 10]. The key step is to view each triangle of spins 
as the seed of a chiral topological phase, a puddle encircled 
by an edge state, as illustrated in Fig. 1. The natural candi- 
date for the chiral topological phase filling the puddle is the 
bosonic v = 1/2 Laughlin state [ ], which is known to have 
the SU(2) symmetry required by our construction and is also 
the state that was envisioned by Kalmeyer and Laughlin [ ]. 
Forming a lattice out of the elementary triangles, we should 
then consider a situation with many individual puddles of this 
topological phase. To see what collective state is formed, we 
have to understand how two puddles are joined. In doing so, 

FIG. 1. (Color online) Sketch of a puddle of topological phase re- 
placing each triangle of three spins. 
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FIG. 2. (Color online) Top panel: Illustration of the behavior of 
the edge states at a corner shared by two puddles. Bottom panel: 
Behavior of two corner-sharing triangular puddles of the topological 
phase. Different colors indicate different chiralities of the topological 
phases. 

it turns out to be useful to consider the specific situation of 
corner-sharing triangles, as realized for instance in a Kagome 
lattice. This situation of edges meeting at the corner shared by 
two triangles is an incarnation of two-channel Kondo physics 
(see e.g. Ref. [33, 34])): If we envision the puddles to be very 
large, they would carry the edge state on each side and the 
corner would look as shown in panel (La) of Fig. 2. The pair 
of edge states on the upper triangle is known [31, 32] to be 
described by the same theory as the right- and left-movers of 
a semi-infinite uniform spin- 1/2 Heisenberg chain, and analo- 
gously for the lower pair of edge states. The spin at the cor- 
ner then appears as the center spin of an infinite chain (pan- 
els (LbJ.c)). It is well known in the context of two-channel 
Kondo physics that the infinite chain will heal if the cen- 
ter spin is coupled to the two semi-infinite chains with equal 
strength [35, 36]. Then, the right- and left-movers will extend 
throughout the entire, infinite system (panel (I.d)). The effect 
on the corner spin is summarized in panel II of Fig. 2, where 
the situation shown in II. a corresponds to La, while Il.d cor- 
responds to I.d. As is evident from Il.d, the corner spin has 
merged the two triangles to form a larger puddle encircled by 
a single edge state, i.e. to form a larger region of the topolog- 
ical phase. 

If we consider all triangles to have the same chirality ( Ja = 
J v ), we can repeat these steps for all corner-sharing triangles 
of the Kagome lattice. The system then forms one macro- 
scopic, extended region of a single topological phase with one 
edge state encircling its boundary, as illustrated in Fig. 3, with 
closed loops encircling the interior hexagons of the Kagome 
lattice. We conclude that choosing uniform chiralities in 



FIG. 3. (Color online) Visualization of the topological phase for uni- 
form choice of chiralities. A collective edge state encircles the whole 
systems, and closed edges encircle each hexagon. 




FIG. 4. (Color online) Real-space (left) and momentum-space (right) 
picture of the staggered (gapless) phase. The colors correspond to 
the three independent directions of extended edge states. The inner 
hexagon, enclosed by a solid line, shows the first Brillouin zone. 

our model, we obtain a direct realization of the Kalmeyer- 
Laughlin proposal for a chiral topological spin liquid phase. 

It is now natural to ask whether other phases can be pre- 
dicted from the same perspective when choosing different pat- 
terns of chiralities. It turns out that this is indeed the case for 
staggered chiralities, i.e. when the down-pointing and the up- 
pointing triangles of the Kagome lattice are assigned opposite 
chiralities ( Ja = — Jy). For this pattern, we find the sought- 
after non-Fermi-liquid state with gapless spin excitations on 
lines in momentum space. 

To see this, we again start from the case of two adjacent tri- 
angles, which are now filled with topological phases of oppo- 
site chirality (Fig. 2 in. a). To see the effect on the edge states, 
consider changing the chirality of the upper triangle by ex- 
changing its two top spins, that is twisting the upper triangle. 
The resulting pattern of edge states after healing is depicted in 
panel ffl.d; note that no backscattering occurs. Again repeat- 
ing this for all triangles of the Kagome lattice results in edge 
states extended throughout the entire system, as shown in the 
left panel of Fig. 4. The pattern of edge states follows three 
sets of parallel lines rotated with respect to each other by 120 
degree rotations. 

To infer the shape of the resulting Bose surface of gapless 
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excitations (corresponding to the extended edge states), we 
note that a single extended edge state, say in the x direction, 
is gapless at a single point in momentum space, k x = 0. Since 
all edges are decoupled, the gapless points corresponding to a 
stack of parallel edge states form a continuous line in momen- 
tum space perpendicular to the direction of k x , i.e. at k x — 
and arbitrary k y . This describes a line of gapless modes in 
momentum space. The full Bose surface for the Kagome lat- 
tice is then obtained when combining three such stacks of in- 
dependent, gapless edge states, each rotated by 120 degrees 
with respect to each other. The resulting Bose surface, shown 
in the right panel of Fig. 4, is therefore comprised of three 
gapless lines [37]. 

Majorana model- Before discussing the technically challeng- 
ing case of SU(2) spins, we illustrate the validity of the above 
network model picture in the situation where the spins are re- 
placed by Majorana fermion zero modes. In particular, we can 
form a term analogous to the spin chirality (1) which defines 
a notion of chirality for a triangle formed by sites i, j, k: 

Xijk = i(iaj + lilk + 7fc7i)- (3) 

Here, ji = 7} denote Majorana operators (7? = 1, {7^ 7^} = 
5ij) on the sites of the Kagome lattice and the choice of sign 
for Jijk again sets the chirality for a triangle. This is a non- 
interacting system which can be solved exactly. 

Similarly to the situation of spins, we can view each tri- 
angle as the seed of a topological phase and apply arguments 
analogous to those discussed in the context of Fig. 2 [36]. Val- 
idating our discussion in the context of Figs. 3 and 4, we find 
that if the chiralities are chosen the same on all triangles, a 
gapped phase emerges (for details, refer to the supplemental 
material [ 38]). Since it is a non-interacting fermion model, the 
Chern number C [ ] is easily calculated [40, 41] and we find 
that C = ±1 depending on the overall chirality of the Hamil- 
tonian. If, on the other hand, we choose staggered chiralities, 
we find a gapless system with a Fermi surface . The shape of 
the Fermi surface precisely matches the discussion above. [ ] 
- Related models have previously been considered [ ] on the 
triangular lattice [43, 44] and Kagome lattice [45, 46] and 
both a topological state and a Fermi surface state have been 
obtained depending on whether the chiralities on all triangles 
are the same or staggered. 

Numerical methods.- Returning to the original spin model (2), 
which is not exactly solvable, the nature of the phases is much 
harder to establish and one has to resort to large-scale nu- 
merical simulations. We briefly summarize our numerical ap- 
proach for (a) the topological phase and (b) the gapless phase. 

(a) Topological phase: To identify the topological nature 
of the phase, we rely on observing the ground-state degen- 
eracy and the presence of edge states. Using exact diago- 
nalization (ED) and the density-matrix renormalization group 
(DMRG) [47-50], we show that there is a gap to all excita- 
tions and that when placed on a torus, there are two low-lying 
states in the S z = sector, which are well-separated from the 
rest of the spectrum. We furthermore show that on a thin, long 
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FIG. 5. (Color online) In all panels, N is the total number of lattice 
sites. Top left panel: Spin gap obtained by DMRG calculations on the 
cylinder. Top right panel: Exact diagonalization gaps of the uniform 
case on both the torus and the cylinder. Bottom panel: Finite-size 
entropy for cylinders. This plot analyzes the entropy at the center of 
the system (averaged over a number of sites to remove oscillations) 
as a function of the total system size. 



strip, the system displays a gapless edge state (see Fig. 3 and 
the supplemental material [38]) which disappears when the 
long boundaries are connected to form a cylinder. All of these 
observations provide evidence for with a bosonic v = 1/2 
Laughlin state. 

(b) Gapless phase: A characteristic of gapless phases with 
a Bose surface is that their quasi-one-dimensional precursors 
are critical with the number of gapless degrees of freedom 
growing with the width of the system [13, 14]. We here focus 
on the particular case of a two-leg ladder, in which two gapless 
modes occur. We therefore expect to observe a central charge 
of c = 2 (see supplemental material [38]). Using DMRG, 
we show that the system is indeed gapless (with the finite-size 
spin gap, A s = E(S Z = 1) + E{S Z = -1) - 2E{S Z = 
0) vanishing as TV -1 , where N is the total number of lattice 
sites) and confirm the expected central charge using fits to the 
entanglement entropy (see supplemental material [38]). 

The DMRG method scales exponentially in the width of the 
system W, but only polynomially or linearly in the length L of 
the system for critical and gapped systems, respectively. This 
allows us to approach thin, long systems, i.e. W small and 
L — > 00. It is based on a variational ansatz that can be system- 
atically improved by increasing the number M of states kept, 
where the computational cost is known to grow as 0(M 3 ). 
We fix W = 2 and use up to M = 3600 states. 
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FIG. 6. (Color online) Top panel: Spin gap for the uniform phase 
on a thin, long strip. Bottom panels: Entanglement entropy for the 
uniform system on a strip. Bottom left: Entropy at the center of the 
system on a semi-logarithmic scale. Bottom right: Central charge 
extracted from a direct fit to the entropy for a given system size. 
Details of the fits can be found in the supplemental material [38]. 



Numerical results: (a) Topological phase.- We start by dis- 
cussing our results for the topological phase obtained for 
uniform chiralities. Our results for the finite-size spin gap 
A S {N), where N is the total number of sites, are shown in the 
upper left panel of Fig. 5. It clearly approaches a finite value 
of A s (oo) « 0.05. To further corroborate that the model has 
a finite correlation length, we can study the scaling of the en- 
tanglement entropy with system size (lower panel of Fig. 5). 
We show clearly that the entanglement entropy at the center 
of the system becomes independent of the total system size 
even for relatively small systems, which is a clear signature 
of a gapped state. This is also consistent with the behavior of 
correlation functions, which we find to decay exponentially 
(see supplemental material [38]). 

To identify the topological nature of this gapped phase, we 
study the dependence of the low-lying spectrum on the topol- 
ogy. We expect a two-fold degeneracy on the torus which is 
split by an exponentially small amount, and no degeneracy 
on the cylinder (i.e. a system with periodic boundary condi- 
tions in the transverse direction, but open boundary conditions 
in the longitudinal direction). The top right panel of Fig. 5 
clearly reproduces this behavior: on the cylinder, we observe 
excited states close to the extrapolated spin gap A s (oo). On 
the torus, on the other hand, we find an additional low-energy 
state in the same S z sector very close to the ground state. 

A final piece of evidence is obtained by placing the system 



on a fully open geometry, i.e. a strip. For this geometry, we 
expect a pair of extended edge states in the long direction to 
give rise to a gapless mode with central charge c — 1. This 
is nicely confirmed by our results in Fig. 6. In the top panel, 
we show the scaling of the spin gap with system size. While 
there is a strong even-odd effect, we clearly observe a scaling 
of A S (N) ~ a/N + b/N 2 (where a, b are fit parameters) as 
expected for a critical system. To measure the central charge, 
we perform fits to the entanglement entropy, as shown in the 
bottom panels of Fig. 6. We find that by either fitting the en- 
tropy at the center of the system as a function of system size, 
or by fitting the entropy as a function of block size for fixed 
system size, we obtain reasonable agreement with c = 1. 
(b) Gapless phase.- Turning to the gapless spin liquid phase 
for staggered chiralities, we again consider the scaling of the 
finite-size spin gap, shown in the top left panel of Fig. 5. We 
find good agreement with A S (N) ~ iV _1 , strongly indicative 
of a critical system. Determining the central charge of the crit- 
ical theory for this gapless system, we find a fit to the entropy 
at the center of the system to be most reliable. Our data, which 
is shown in the lower panel of Fig. 5, shows good agreement 
with the expected value c = 2 of the central charge. 
Conclusions - Summarizing our results, we have obtained 
strong numerical evidence for (i) a gapped topological spin 
liquid and (ii) a gapless spin liquid with gapless spin excita- 
tions on lines in momentum space. These phases are numeri- 
cally observed in a relatively simple spin- 1/2 model of three- 
spin interactions around the triangles of the Kagome lattice. 
The key concept to understand the emergence of these phases 
is a network model of edge states. For the case of the gap- 
less spin liquid, the geometry of the gapless lines in momen- 
tum space follows directly from the aforementioned network 
model picture. The same reasoning, which leads from the 
microscopic network model to the emergent nature of these 
phases, can also be used and thereby confirmed in a differ- 
ent, exactly solvable model where the spins are replaced by 
Majorana fermion zero modes on the same lattice. 

One aspect of our work that we plan to address in the future 
is the question of interactions between the gapless excitations 
on the lines in momentum space in the spirit of generaliza- 
tions of Landau's Fermi liquid theory [ I ]. We also plan to 
explore the phase diagram when additional terms, such as a 
Heisenberg term on the bonds, are added to our model. 

Future work should explore whether decoupling the 
Heisenberg antiferromagnet on the Kagome lattice using the 
well-known relation xfjk ~ (&i + &j + $k) 2 [21] can lead to 
new insights into the much debated nature of the ground state 
of this model. Another interesting future direction to pursue is 
to make a connection to orbital currents induced by finite spin 
chiralities, as argued in Ref. [52]. 
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SUPPLEMENTARY MATERIAL 

This supplementary material is divided into three sections: 
in Section A., we discuss in more detail how the gapless lines 
in momentum space are obtained from the network model of 
edge states for the model with staggered chiralities (gapless 
spin liquid). In Section B., we explain in detail how the quasi- 
one-dimensional precursor systems of the two-dimensional 
gapped and gapless phases can be understood in terms of the 
network picture. In Section C, we focus on the special case 
of the (exactly solvable) Majorana hopping problem, which is 
obtained when the spins of the original model are replaced by 
Majorana zero modes, and we give some detail on its diago- 
nalization. In particular, we calculate the shape of the Fermi 
surface explicitly for this model. In Section D. we briefly dis- 
cuss entanglement entropy and central charge, whereas in Sec- 
tion E. we discuss the exponential decay of spin correlation 
functions in the gapped, topological phase. 



A. Gapless lines in momentum space 

In Figure 1 0, we show in more detail how to infer the loca- 
tion where gapless states exist in momentum space (i.e, Fermi 
surface in the case of fermions) from the network model of 
edge states in real space. The top left panel of Fig. 10 shows 
a stack (in the y direction) of uncoupled edge states directed 
in the x direction. The resulting picture in momentum space 
is shown in the lower left panel of the figure: An edge state 
in the (say) x direction is gapless at one point in momentum 
space, i.e. when the momentum k x vanishes, k x = 0. Since 
the different edges are decoupled, the gapless points of the 
stack of parallel edge states form a continuous line in momen- 
tum space perpendicular to the direction of the edges, i.e. at 
k x = and arbitrary k y G [—%,%). This describes a line of 
gapless modes in momentum space. 

The full real-space network on the Kagome lattice, which 
is shown in Fig. 4 of the main paper, is recovered when three 
such stacks of gapless edges, rotated by 120 degrees with 
respect to each other, are combined. In the center panel of 
Fig. 10, we show the situation for two such stacks, i.e. a sec- 
ond stack of such edges, rotated by 120 degrees with respect 
to the edges in x direction, has been added. This additional 
stack of edges will give rise to an additional line of gapless 
states in momentum space, which is rotated by 120 degrees 
with respect to the first one as well, as shown in the lower cen- 
ter panel of Fig. 10. In the right-most panel, all three stacks 
are combined, giving rise to three lines of gapless states in 
momentum space rotated with respect to each other by 120 
degrees. 



B. Quasi-one-dimensional precursors of the two-dimensional 
phases 

In this section, we provide additional information on 
the quasi-one-dimensional precursors of the two-dimensional 
phases. Since the DMRG method scales exponentially in the 
width of the system while scaling only polynomially in its 
length, we perform our numerical simulations on quasi-one- 
dimensional systems. 

Figure 1 1 shows a quasi-one-dimensional version of the 
networks discussed in the main text. The top panel refers 
to the case where the corresponding two-dimensional system, 
depicted in Fig. 3 of the main text, is the gapped 2D phase. 
The bottom panel refers to the case where the correspond- 
ing two-dimensional system, depicted in Fig. 4 of the main 
text, is the gapless 2D phase. In the figure, we show a sys- 
tem with open boundary conditions in the transverse direc- 
tion, and with a width of W = 2 and a total number of sites 
N = 40. Periodic boundary conditions in the transverse di- 
rection are obtained by identifying the top-most spin of the 
top row of (up-pointing) triangles with the bottom-most spin 
of the bottom row of (down-pointing) triangles. In this case, 
the topology is that of a surface of a cylinder or torus, for open 
or periodic boundary conditions in the longitudinal direction, 
respectively. 

The upper panel of Fig. 1 1 shows the uniform case, which 
is predicted to be in a topological gapped chiral spin liquid 
phase. In the bulk, this phase has no extended edges; in- 
stead, all edges combine to form localized states encircling 
the hexagons. These are shown as green dashed lines in the 
figure. For the open boundary conditions shown in the fig- 
ure, extended edge states form along the boundaries of the 
system in the long (horizontal) direction. These edge states 
are shown as red dashed and blue dotted lines in the figure. 
From this picture, we predict that the central charge for this 
quasi-one-dimensional system is c = 1 (corresponding to the 
edge of a single bosonic v = 1/2 Laughlin state described 
by an SU(2)i CFT), which is confirmed by the numerical cal- 
culation of the entanglement entropy, which is detailed in the 
main paper. 

The lower panel of Fig. 1 1, on the other hand, shows the 
staggered case, which is predicted to be a gapless spin liq- 
uid phase. For the quasi-one-dimensional system shown in 
the figure, one set of extended edges forms along the horizon- 
tal lines of the Kagome lattice, shown here as solid red line. 
These are right-moving. While in the truly two-dimensional 
system, two additional, independent stacks of extended edges 
would appear on the two other sets of lines which are rotated 
by 120 and 240 degrees, the quasi- Id geometry forces these 
edges to follow the zig-zag pattern shown as dotted blue lines 
in the figure. Those edges are left-moving. Counting the total 
number of edges shows that there are two left-moving and two 
right-moving edges leading a the central charge of c = 2. 

Figure 7 shows how the results for the quasi-one- 
dimensional systems discussed in the previous paragraph in- 
dicate the existence of extended lines of gapless modes in mo- 
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FIG. 7. (Color online) Gapless lines (Fermi surface) of the two- 
dimensional model, cf. Fig. 4 of the main text. In this figure, we ad- 
ditionally sketch a set of allowed momenta (dashed blue lines) which 
are obtained if the system is placed on a quasi-one-dimensional ge- 
ometry. 

mentum space in the two-dimensional gapless phase. Making 
the system finite in the transverse direction leads to discretized 
momenta in this direction, where the number of such allowed 
momenta grows linearly in the width of the system. In Fig- 
ure 7, these allowed momenta are indicated as blue dotted 
horizontal lines. The system is now gapless only at a finite 
number of points in momentum space which are given by the 
intersection of the allowed momenta with the gapless lines of 
the fully two-dimensional system. For such a gapless phase, 
the central charge measured in a quasi-one-dimensional ge- 
ometry grows linearly in the number of such intersections. 

C. Non-interacting Majorana model 



FIG. 8. (Color online) Intersection of the central band E2(k x , k y ) 
of the staggered model with the E — plane (light blue). The lines 
where the E = plane crosses the band represent the gapless lines; 
areas with E2(k x , k y ) < correspond to filled regions of momen- 
tum space. 

of states in this band coincides with the number of hexagons 
in the system, we identify these zero-energy states as the non- 
interacting edge states encircling the hexagons of the Kagome 
lattice, which are shown in Fig. 3 of the main text. We calcu- 
late the Chern number using the prescription of Refs. [40, 41] 
and find that the Chern number C of the top and bottom band 
is C = ±1, i.e. that the model is in a topological phase. 

Our results for the staggered phase are summarized in Fig- 
ures 12 and 8. Fig. 12 shows all three bands. The top band has 
Es(k x , k y ) > 0, while the bottom band has E%(k x , k y ) < 0. 
All three bands are gapless and touch at k x = k y = 0. In 
Fig. 8, we show the intersection of the central band E% with 
the E — plane. This clearly shows the three straight lines 
where this band becomes gapless, as well as the regions of 
momentum space of filled and empty states. The agreement 
between this and the picture developed in the main text con- 
firms that our reasoning, which is based purely on considering 
a network model of edge states, is correct. 



In the following section, we will discuss some details of the 
band structure obtained for the Majorana hopping Hamilto- 
nian (cf. Eqn. (2)) 

H = Xijk ± X! Xijk (4) 
a v 

Xijk = AlH] + Ijlk + Ikli), (5) 

where the case with the + (— ) sign corresponds to the uniform 
(staggered) phase. In both cases, the Hamiltonian is quadratic 
in the Majorana fermion operators and can therefore be diag- 
onalized straightforwardly. Since the Kagome lattice is ob- 
tained as a triangular lattice of three-site unit cells, we obtain 
three energy bands E a (k x , k y ), a = 1, 2, 3. 

In the uniform phase, we observe that all three bands are 
separated by a gap from each other, and the central band E% 
is dispersionless, i.e. Ei(k x , k y ) — 0. Noting that the number 



D. Entanglement entropy and the central charge 

From the ground states obtained by DMRG, the entangle- 
ment entropy for a contiguous block of n sites at the end of 
an open system of N total sites can be extracted without ad- 
ditional cost. This allows us to extract the central charge by 
performing a fit to the well-known result [54, 55] for the en- 
tanglement entropy of a contiguous block of sites in an open 
Id system, 

S(n)=5Q + |log^ S m^V (6) 
which for n = NJ 2 reduces to 

S(n = N/2) = S + C - log (—). (7) 
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FIG. 9. (Color online) Correlation function in the topological phase 
for a cylindrical topology with W = 2. The correlation function 
shown is measured along one of the two horizontal chains that com- 
prise the system, and i denotes the distance from the end of the sys- 
tem along the chain. 



Using Eqn. (6), the central charge can be obtained from the 
ground state for a given system size N by fitting to the entan- 
glement entropy for various subsystem sizes n. In most cases, 
this is the most efficient way to extract the central charge. 
In some cases, however, it is more reliable to use Eqn. (7), 
which requires the ground state for several different total sys- 
tem sizes N. 

E. Exponential decay of correlation functions in the topological 

phase 

In Fig. 9, we show a correlation function in the topological 
phase. We have measured (SqS~), where all sites are chosen 
to be on one of the horizontal chains that comprise this system, 
and i measures the distance along this chain. We find very 
good agreement with an exponential decay, i.e. 

|(S+S-)|~expH/£) (8) 

with £ w 0.44. This exponential decay is consistent with a 
gapped, topological phase. 
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FIG. 10. (Color online) This figure illustrates the construction of the gapless lines in momentum space (lower panels) from the real-space 
picture (upper panels). Each of the three sets of parallel lines forming the Kagome lattice gives rise to a line of gapless points in momentum 
space. The horizontal lines give rise to gapless states at k x — (left panel). The two other lines are rotated by 120 and 240 degrees, 
respectively. 
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FIG. 11. (Color online) Two-leg Kagome systems with illustration of the choice of chiralities for the uniform and staggered model, in quasi- 
one-dimensional stip geometry. The red dashed lines indicate right-moving extended states, blue dashed line indicate left-moving extended 
states, and the green dashed lines encircling the hexagon show the localized states in the topological phase. 




